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ABSTRACT 

Using the Parkes radio telescope we have carried out deep observations of eleven 
unassociated gamma-ray sources. Periodicity searches of these data have discovered 
two millisecond pulsars, PSR J1103-5403 (1FGL J1103.9-5355) and PSR J2241-5236 
(1FGL J2241.9-5236), and a long period pulsar, PSR J1604-44 (1FGL J1604.7-4443). 
5^ , In addition we searched for but did not detect any radio pulsations from six gamma- 

ray pulsars discovered by the Fermi satellite to a level of ~ 0.04 mJy (for pulsars with 
a 10% duty cycle). 

Timing of the millisecond pulsar PSR J1103-5403 has shown that its position is 
9' from the centroid of the gamma-ray source. Since these observations were carried 
out, independent evidence has shown that 1FGL J1103. 9-5355 is associated with the 
flat spectrum radio source PKS 1101-536. It appears certain that the pulsar is not 
associated with the gamma-ray source, despite the seemingly low probability of a 
chance detection of a radio millisecond pulsar. We consider that PSR J1604-44 is a 
chance discovery of a weak, long period pulsar and is unlikely to be associated with 
1FGL J1604.7-4443. PSR J2241-5236 has a spin period of 2.2 ms and orbits a very 
low mass companion with a 3.5 hour orbital period. The relatively high flux density 
and low dispersion measure of PSR J2241-5236 makes it an excellent candidate for 
high precision timing experiments. The gamma-rays of 1FGL J2241. 9-5236 have a 
spectrum that is well modelled by a powerlaw with exponential cutoff, and phase- 
binning with the radio ephemeris results in a multi-peaked gamma-ray pulse profile. 
Observations with Chandra have identified a coincident X-ray source within 0.1"of the 
position of the pulsar obtained by radio timing. 
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1 INTRODUCTION 



The Large Area Telescope (LAT, lAtwood et all 120091 1 
aboard the Fermi Gamma-ray Space Telescope has detected 
more than six hundred unassociated gamma-ray sources 
amongst the ~ 15 00 sources in the f irst year point source 
catalogue (1FGL. IXbdo et all l2010ah . Of the gamma-ray 
point sources which are associated with known astronomi- 
cal objects, almost all fall into one of two group s: radio-loud 
nuclei of distant galaxies and energetic pulsars (|Abdo et al.l 
l2010bl R>. In the last 15 years there has been great inter- 
est in observing unassociated gamma- ray sources with the 
aim of detecting radio pulsars (e.g. Roberts et al. 2002; 
Keith et all 120081 ; ICamilo et all 120091; lAbdo et all l2010dl ; 



Rav et alfeOlOf ). Prior to the launch of Fermi, the positional 



error-circles of gamma-ray sources have generally been larger 
than the primary beam of large radio telescopes (at the typi- 
cal observing frequency of 1400 MHz). Because of this, these 
searches have required many pointings to cover each source, 
limiting the sensitivity of the searches and the number of 
sources that can be covered. 

The increased sensitivity and spatial resolution of the 
LAT means that the error circles of most sources are com- 
parable to the typical radio beam size, allowing for much 
deeper radio searches with a single pointing per source. The 
radio pulsar community is conducting a coordinated effort 
to observe the unassociated gamma-ray sources in the radio 
band. There are two main benefits for identifying gamma- 
ray sources as radio pulsars. Firstly, we can remove these 
objects from the unassociated gamma-ray list and therefore 
hope to better understand of the nature of the remaining 
unassociated sources. Secondly, in most cases, the radio de- 
tections are essential for studying the gamma-ray pulsars 
themselves. Radio observations of a pulsar can give a pre- 
cise ephemeris for phase binning the gamma-ray photons 
to form a gamma-ray pulse profile. The phase-resolved in- 
tensity and spectra of gamma-ray pulsars are providing a 
great deal of insight into the pulsar magnetosphere and al- 
lowing for distinction between the various models of gamma- 
ray emission, especially w hen compared with radio profiles 
jRomani 

Additionally, the LAT sensitivity has allowed for the 
discovery of 24 pu l sars by their gamma-ray e mission alone 
l|Abdo et al]|2009al ; ISaz Parkinson et al.ll2010f) . Understand- 
ing the fraction of gamma-ray pulsars with detectable radio 
emission has implications for radio and gamma-ray emis- 
sion models (|Ravi et al.ll2010l V Deep observations of these 
gamma-ray pulsars are therefore required to confirm or rule 
out pulsed radio emission. 

The pulsar population is typically categorised based on 
the observed spin period, P, and its derivative, P. There 
are two classes of pulsars for which the rate of rotational 
energy loss, E oc PP~ 3 , is high enough to enable gamma- 
ray emission. The most energetic are 'young' pulsars, such 
as the Vela pulsar (PSR J0835-4510), which typically have 
spin periods between 30 and 300 ms, and period derivatives 
above 10 -15 ss _1 . In contrast, millisecond pulsars (MSPs) 
such as PSR J0218+4232, have very small period derivatives 
but high E due to their fast rotation periods, often below 
10 ms. Whilst E appears to give a good measure for potential 
for gamma-ray emission, sensitivity is also limited by the 
pulsar distance. 



This paper presents a search with the Parkes radio tele- 
scope for young pulsars and MSPs, targeting unassociated 
gamma-ray sources and gamma-ray pulsars discovered by 
the Fermi telescope. 



2 SOURCE SELECTION 

In March of 2009 we selected sources from a preliminary 
list of Fermi point sources for which no formal association 
with a known astronomical object had been made. To best 
suit observation with the Parkes radio telescope, we selected 
sources with a declination south of —40° and an error circle 
with diameter smaller than the 14 arcminute beam width 
of Parkes at 1400 MHz. These criteria allowed us to select 
11 sources, as shown in Table [1] As the Fermi LAT all- 
sky survey continues, gamma-ray photon statistics increase 
and localisation improves. Therefore the positions observed 
in this search are offset somewhat from the more accurate 
position as published the 1FGL catalogue and these offsets 
are shown in the table. Additionally, we chose to observe 
six southern gamma-ray pulsars that had not been deeply 
searched using radio telescopes, with the aim of detecting 
any radio pulsations. The observations of these pulsars are 
detailed in Table H 



3 OBSERVATION AND DATA ANALYSIS 

Each of the selected sources was observed using the Parkes 
64 m antenna and the centre beam of the 21-cm 'Multibeam' 
receiver, which has a noise equivalent flux density of ~40 Jy. 
The BPSR digital interbank system was used which has an 
observing bandwidth of 340 MHz centred at 1352 MHz and 
split into 870 channels. Each channel is sampled with 2-bits 
of precision every 64 [is. In order that we might place strin- 
gent limits on the detectability of radio pulses from these 
objects, we observed for 4.8 hours per target. We determine 
the minimum detectable flux density according to the stan- 
dard interpretation of the radiometer equation, 



Smin — 



fj(T sys + T s ky) / Wf vac 



(1) 



For this paper we assume: system temperature, T sys = 23 K; 
sky temperature T s k y = 2.5 K; gain, G = 0.735; bandwidth 
B — 340 MHz; observing time r bs = 17200 s; effective pulse 
duty cycle, Wf ra c = 0.1; and the signal-to-noise threshold, 
a = 8. For a pulsar close to the centre of the beam, this 
gives a detection threshold of ~ 0.03 mjy. A detailed dis- 
cu ssion of the s ensiti vity of the BPSR system can be found 
in lKeith et all (|201Ch . 

For the known gamma-ray pulsars we folded the data 
using a rotational ephem eris provided by the Fermi collab- 
oration (|Rav et al.l |2010| ) , searching in dispersion measure 
(DM) from to 1000 cm -3 pc. The data taken on unasso- 
ciated 1FGL sources were searched for periodic signals us - 
ing the hitrun pipeline as described in lKeith et al.l l|2010h . 
Again, DMs up to 1000 cm -3 pc were searched. We have 
searched the data attempting to correct for the period drifts 
due to acceleration of the pulsar if in a binary system. The 
acceleration searching was performed using presto, using 
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Table 1. List of the Fermi sources targeted in this search. Columns are the 1FGL catalogue name, the coordinates of the observation 
(J2000), the offset from the 1FGL position, the date of observation. The rotation period and dispersion measure are listed for the three 
pulsars discovered in this search + . 
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RA 
(h 


m 
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' ") 
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(') 


Obs Date 


Period 
(s) 


DM 

(cm -3 pc) 


J0823. 


.3-4248 


08: 


:23: 


:26. 
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-42 


:50 


:38 


,4 


2. 


2 


2009-08-03 






J1018. 


.6-5856 


10: 


: 18: 


:30. 


2 


-58 


:57 


:47 


.5 


2. 


1 


2009-04-14 






J1045. 


2-5942 


10: 


:45: 


:11. 


1 


-59: 


:43 


: 10 


,2 


0. 


9 


2009-04-15 






J1103. 


9-5355 


11: 


:03: 


:53. 


8 


-53: 


:55 


:58 


,8 


0. 


3 


2009-08-01 


0.0033927 
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J1514. 


.1-4945+ 
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:44 


:30 


,1 


2. 


3 


2009-04-13 






J1536. 


.5-4949 
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:36: 


:26. 


4 


-49 


:47 


:48 


.5 


1. 


9 


2009-04-15 






J1604. 


.7-4443 
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:32. 


9 


-44 


:42 


:43 


.2 


2. 





2009-08-01 


1 . 3892 
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J1620. 


8-4928C 
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:05. 
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-49: 


:30 


:32 


,4 


3. 





2009-08-03 
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16: 


:35: 
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5 


-47: 


:28 


:55 


,2 


13. 





2009-04-13 






J1640. 


. 8-4634C 


16: 


:39: 


:59. 





-46 


:44 


:24 


.0 


13. 





2009-08-02 






J2241. 


.9-5236 


22: 


:41: 


:55. 
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-52 


:36 


:39 


.6 


0. 


3 


2009-08-01 


0.0021866 


11 



+ A pulsar was discovered in 1FGL J1514. 1—4945 in another search at Parkes, as discussed in Kerr et al. (in prep). 



the algorithm described in iRansom et all (|2002l ). The fre- 
quency drift range searched was ±150 Fourier bins, corre- 
sponding to an acceleration range of ±0.3 ms~ for a pulsar 
with a rotational frequency of 500 H j]]. Additionally we have 
done the same processing on the first 2100 s of each obser- 
vation, giving an acceleration range of ±20.4 ms~ 2 for the 
same pulsar but a reduced sensitivity of 0.1 mjy. Further 
processing to remove the effects of binary motion will likely 
be done in the future. From the candidates presented by the 
software, three were selected for re-observation at Parkes, 
and have been confirmed to be radio pulsarfl 



3.1 Simulations 

Naively, one might assume that there is little chance of de- 
tecting rare objects such as MSPs in so few pointings, unless 
the MSP is indeed the source of the gamma rays. In order to 
determine the probability that any of our detections could be 
unrelated to the gamma-ray sources, we chose to simulate 
the Galactic pulsar populatio n with the PSRPOP software, 
using the method described in iKeith et al l (|201Ch . based on 
the model of lLorimer et alj (|2006t ). For the non-MSP pop- 
ulation our model has a log-normal spin period distribu- 
tion (with mean and variance in log-space of 2.71 and 0.34 
respectively), a power-law luminosity distribution (with in- 
dex —0.59, cutting off below 0.1 mjy kpc 2 ), an exponential 
distribution for the height above the Galactic plane (with 
scale 0.33 kpc) a n d usi ng the radial distribution given by 
lYusifov fc Kuciikl (|2004l ). For the MSPs, we used a similar 
model except that, to better match the known MSP popu- 
lation, the scale height was increased to 0.5 kpc and the pe- 
riod distribution was empirically matched to the population 
of known MSPs. The total number of pulsars simulated was 
limited such that the number of pulsars d iscovered in a sim- 
ulate d Parkes Multibeam Pulsar Survey ^Manchester et all 
l200ll ) matched the real discovery count. 

We then simulated the 11, 4.8-hour pointings of our 



1 The acceleration range goes as l/i/ 

2 A fourth pulsar, detectable in our observation of 
1FGL J1514. 1—4945, was discovered in the search described by 
Kerr et al. (in prep). 



search that did not point at known gamma-ray pulsars, and 
recorded the number of pulsars discovered. This was re- 
peated 1000 times with new simulations of the Galactic pul- 
sar population each time, based on different random seeds. 
The mean detection rates were 1.3 non-MSPs and 0.05 MSPs 
in 11 pointings. Therefore it is likely that we could discover 
at least one non-MSP in our searches, and the discovery of 
an unassociated MSP is certainly not ruled out. 



4 RESULTS AND DISCUSSION 

None of the six pulsars discovered by their pulsed gamma- 
rays, as listed in Tabled was detected in our radio searches. 
The discovery of these pulsars and searches for the ir ra- 
dio pulsations are detail ed in lAbdo et all (|2009al ) and 
ISaz Parkinson et al.l (|2010h . Analysis of the other 11 obser- 
vations has yielded the discovery of three radio pulsars, and 
the detection of the MSP J1514-49 (discussed in Kerr et al., 
in prep). Table [1] presents the spin period and DM of the 
pulsars discovered. PSRs J1103-5403 and J2241-5236 have 
spin periods less than 5 ms and therefore belong to the class 
of MSPs. Although we initially expected to discover mainly 
young, energetic pulsars in this search, we are finding that 
many MSPs are being detected in gamma rays (|Abdo et al.l 
1200901 ). PSR J1604-44 has a spin period of 1.4 s, typical 
of old, but not recycled, pulsars which are not commonly 
associated with gamma-ray emission. In this section we dis- 
cuss the non-detections, and present relevant details of the 
gamma-ray and radio analysis of PSRs J1103-5403, J1604- 
44 and J2241-5236. 



4.1 Non-detections 

No radio pulsars were detected in observations of seven of 
the unassociated point sources and all six of the gamma-ray 
pulsars. 1FGL J1636.4-4737c and J1640.8-4634c are con- 
fused with diffuse Galactic emission and are now considered 
to lie a considerable distance from the observed location, 
therefore our radio observations are not constraining. For 
the 11 other sources we are confident that any radio pulsa- 
tions associated with the source must be either be below our 
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Table 2. The six gamma-ray selected pulsars that were targeted 
for radio observations as part of this work. 
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2009-04-14 


J1809- 
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09 
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-23 


32 
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2009-04-15 



minimum detectable flux density, hidden by effects of the 
interstellar medium (scattering or scintillation), or by the 
effects of binary motion. For the sources for which gamma- 
ray pulses have not been detected, our Fourier based search 
method is only maximally sensitive to solitary pulsars, or 
those with orbital periods of greater than a few days. Al- 
though our search has used some acceleration compensation, 
for pulsars with large changes in the apparent pulse period 
our sensitivity will be considerably reduced. 



4.2 PSR J1103-5403 and 1FGL J1103. 9-5355 

4.2.1 PSR J 11 03-5403 

PSR J1103-5403 is a solitary pulsar with spin period of 
3.4 ms and DM 104 cm -3 pc. Timing parameters are given 
in TableO These parameters indicate that PSR J1103-5403 
is a MSP at a distance of 2.5 kpc, u sing the electron den- 
sity model of ICordes fc Laziol (|2002T ). The average profile 
appears relatively simple, with a small amount of linear and 
circular polarisation, as shown in Figure [1] 

Intriguingiy, since its selection for this survey, the 
gamma-ray source 1FGL J1103.9— 5355 has been associ- 
ated with the flat s pectrum radio source PKS 1101-536 
l|Mahonv et alj|2010h . believed to be an active galactic nu- 
cleus (AGN). Since both MSPs and AGNs are strong candi- 
dates for gamma-ray emission, further analysis is required to 
determine which is the true source of 1FGL J1103.9— 5355. 
Since the discovery, our timing observations constrain the 
position of the pulsar to Ilh03m33s — 54°03'43", well out- 
side the 2.4' error radius (95% confidence) of 1FGL J1103.9- 
5355. Nevertheless, as described below, we have performed 
a careful analysis of the gamma-ray data to confirm that the 
sources are distinct. 



4-2.2 Gamma-ray data selection 

We used 21 months of Fermi LAT data taken in sky-scanning 
survey mode between 4 August 2008 - 4 May 2010 (MET0 
239557417- 294680617). This extends the data set used for 
the 1FGL catalogue by an additional ten months. As in the 
1FGL catalogue, we used only pass 6v3 diffuse class photons. 
We selected events within a 10° radius of 1FGL J1103.9- 
5355 in the energy range (100 MeV-300 GeV). The software 
used for the analysis was the LAT Science Tools v9rl5p2, 



Table 3. Astrometric and rotational parameters from the timing 
of PSR J1103— 5403. Values in parentheses are the nominal la 
uncertainties in the least significant digits quoted. 



PSR J1103-5403 

Right ascension, a 11:03:33.294(3) 

Declination, S -54:03:43.245(16) 

Dispersion measure, DM (cm _3 pc) 103.915(3) 

Rotation measure, RM (rad m~ 2 ) -103(21) 

Pulse frequency, v (s^ 1 ) 294.749654528(3) 

First derivative of pulse frequency, i> (s — 2 ) — 3.2(8) X 10~ 16 

Epoch of model (MJD) 55300 

Validity of solution (MJD) 55044.1—55459.9 

Rms of residuals (/is) 6.42 

Clock correction procedure TT(TAI) 

Solar system ephemeris model DE405 

Coordinate time standard TCB 



Flux density at 1.4 GHz (Si 40 o, mjy) .... 0.18(4) 

Width at 50% of pulse peak (Wso, /is) ... . 140 

Distance* (kpc) 2.5 

Characteristic age (yr) 1.4 X 10 10 

Surface magnetic field strength (G) 1.1 X 10 8 

Rate of energy loss, E (ergs -1 ) 3.8 X 10 33 

assumes the lCordes fc Laziol (2002) model of electron density. 




-60 -50 -40 -30 -20 -10 10 20 30 40 50 60 
Pulse Phase (deg) 

Figure 1. The 1369 MHz integrated pulse profile of PSR J1103- 
5403 clipped to a region round the pulse peak. Total intensity is 
shown as a thick line, with linear polarisation as a thin solid 
line and circular polarisation as a dotted line. Flux density is 
measured as a fraction of the pulse peak. The polarisation position 
angle variation over the pulse is shown in the upper panel. The 
polarisation position angles shown are not absolute. The inset 
figure shows the total intensity profile over the entire 360°. 



3 Mission Elapsed Time, or seconds since 00:00 UTC on 2001 
January 1 
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Figure 2. The measurements of position of 1FGL J1103.9- 
5355 and nearby Swift sources shown with the radio position of 
PKS 1101-536 and PSR Jl 103-5403. The large circles show the 
95% confidence limit of the position from the 1FGL catalogue, 
'off-pulse', and 'all- phase' analyses as labelled. The two Swift 
sources close to the gamma-ray source are indicated by small cir- 
cles and labelled Source 1 and 2. The cross labelled PKS 1101-536 
shows the position of the known active radio galaxy, and is coin- 
cident with Swift source number 2. The cross labelled PSR J1103- 
5403 marks the position of the pulsar obtained from radio timing. 

available from the Fermi Science Support Centre (FSSC) 
websit^E with the instrument resp onse functions PAS S6-V3. 

The pulsar PSR J1057-5226 l|Abdo et alj|2010eft is lo- 
cated 1.8° from 1FGL J1103. 9-5355 and its large gamma-ray 
flux could therefore have an effect on our spectral and local- 
isation results. In order to remove (or at least reduce) the 
effect of the pulsar, we selected the subset of events falling 
in the off-pulse region of PSR J 1057-5226, according to an 
updated ephemeris for the pulsar derived from LAT data 
covering the entire 21-month period of observations 0. We 
refer to this reduced data set, consisting of ~60% of the 
events, as the 'off-pulse' data set. When the entire data set 
is used, we will refer to this as the 'all-phase' data set. 

4-2.3 Gamma-ray spectrum and localisation 

We performed an unbinned likelihood analysis using the 
LAT Science Tool gtlike, adopting the standard mod- 
els for the diffuse emission components, both Galactic 
(gll_iem_v02 . f it) and isotropic (isotropic_iem_v02.txt). 
All 1FGL sources within 15° of 1FGL Jl 103.9-5355 were 
included in the spectral model. We first assumed a power- 
law spectrum for all sources and then refined the model by 
introducing an exponential cutoff for 12 of them, after de- 
termining, iteratively, by means of a likelihood ratio test, 

4 http : / / fermi .gsfc. nasa.gov /ssc / 

5 The timing model used in this study will be made available on 
the servers of the FSSC. 



that the addition of such a parameter resulted in a signif- 
icant improvement of the model. For all sources within 6° 
of 1FGL J1103. 9-5355 we fitted the normalisation, while for 
PSR J1057-5226 and 1FGL J1117.0-5339 all the spectral 
parameters were allowed to be free. 

The addition of a cutoff in the spectrum of 
1FGL J1103. 9-5355 did not lead to a statistically significant 
increase in the likelihood of the model, and we therefore 
modelled the source with a simple power law. Using the 'all- 
phase' sample, we obtained a best fit value for the spectral 
index of 2.20 ± 0.06, broadly in agreement with what was 
reported in the 1FGL catalogue (2.05 ±0.06). The resulting 
(0.1-300 GeV) photon flux is (6.3 ± 0.8) x 10" 8 cm" 2 s _1 . 
We repeated our analysis using the 'off-pulse' data set and 
obtained a spectral index of 2.20 ± 0.07 and a photon flux 
of (6.2 ±0.8) x 10~ 8 cm _2 s -1 , consistent with the previous 
results. 

Using the LAT Science Tool gtf indsrc we attempted 
to improve on the 1FGL localisation by performing a maxi- 
mum likelihood fit on our longer data set and more refined 
spectral model. We allowed the parameters of sources within 
2° of 1FGL J1103. 9-5355 to vary, while fixing all the other 
parameters to the values obtained from our spectral anal- 
ysis. In order to check the reliability of this technique, we 
used gtf indsrc and our spectral model for the region to 
derive a position for PSR J1057-5226. The resulting best fit 
position is 10h57.96m -52°26.8', with a 95% (statistical) 
uncertainty radius of r95%=0.84', less than 0.18' away from 
its well-determined radio position. 

For 1FGL J1103. 9-5355, using the 'all-phase' data set, 
we obtain a best fit position of llh04.08m — 53°55.6', and 
r95%=1.3'. To check for possible effects from the nearby 
pulsar PSRJ1057-5226 (described in Section l4~2~2"j) we ran 
gtf indsrc on the 'off-pulse' data set and model, obtaining 
a best fit position of llh03.98m -53°54.8', with r95%=2.3', 
consistent with our previous result. 

Using a ~4 ks 0.2-10 keV Swift XRT observation (Ob- 
sID 00031616001), we searched for possible X-ray counter- 
parts of 1F GL J1103.9-5355. Fiv e sources are detected in the 
Swift XRT (jBurrows et alj|2005h field of view, but only two 
are close enough to be consistent with the LAT gamma-ray 
source (labelled 1 and 2 in Figure [5J. Source 1 is located at 
Ilh04m03.37s -53°56'05.50" (with r90%=6.65") and has a 
count rate of 2.63 ± 0.99 x 10 3 cts s _1 , while source 2 is 
located at Ilh03m52.02s -53°56'56.6l" (with r90%=4.76") 
and has a count rate of 13.9 ± 2.1 x 10 3 ctss -1 . Consis- 
tent with the position of Swift X-ray source 2 is a bright ra- 
dio source, PKS 1101-5 36, located at Ilh03m52.221639(13)s 
-53°57'00.69607(20)" l|johnston et al.lll995l ; iMahonv et al.l 
2010). These gamma-ray, X-ray and radio positions are sum- 
marised in Figure [2] 

4-2.4 Variability analysis 

Since gamma-ray pulsars, unlike AGNs, are typically not 
variable on long time-scales we investigate the variability 
characteristics of 1FGL J1103. 9-5355. The 1FGL catalogue 
gives a variability index of V=48.8 for this source, implying 
that it has a proba bility of <1% of being a steady source 
[Abd o et al. I l2010ah . Here, we extend the variability analy- 
sis to our 21-month data set. For this variability analysis, 
we used the 'off-pulse' data set as defined in Section 14.2.21 
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Figure 3. Fermi LAT monthly light curve of 1FGL J1103.9-5355 
in the 0.1-300 GeV energy range. 



to avoid any variability induced due to the fluctuations in 
the spectral fit of the nearby strong pulsar PSR J1057-5226. 
We obtained a monthly light curve for 1FGL J1103. 9-5355 
dividing the data set into 30-day intervals and performing a 
spectral analysis in each time bin. In each time bin, we used 
the same spectral model as described in Section 14.2.31 for 
nearby sources and diffuse components, allowing only their 
normalisation parameters to be free in the fit, and fixing 
their spectral indices to the best fit value obtained from the 
analysis of full time period. We obtained the flux of 1FGL 
J1103. 9-5355 for each time bi n it was detected w ith a test 
statistic, TS, greater than 9 jMattox et al.|[T996l ). and we 
calculated a 95% confidence level upper limit for those time 
bins that the source was not detected. We found that 1FGL 
J1103. 9-5355 is clearly variable as seen from Figure 3; it ex- 
hibits an excess emission over its average flux level early in 
the observation period, while it drops below the detection 
limit in the first two and five of the last 6 months. Since the 
monthly light curve turns into upper limits in some time 
bins, we repeated the spectral analysis splitting the data 
set into 4 segments, each 160 days, sufficiently long to de- 
tect the source in every time bin. We got the following val- 
ues for the photon flux of 1FGL J1103. 9-5355 (in units of 
10" 8 cm" 2 s _1 ): 9.2±1.6, 10.4±1.5, 4.2 ±1.2, and 1.1 ±0.7. 
The corresponding fit to a constant has a x 2 = 45.65 with 
three degrees of freedom, allowing us to reject the hypoth- 
esis of a steady source with a confidence level of > 6a. A 
very conservative estimate of the systematics as 20% of the 
measured flux still confirms the variability of the source at 
a 3er level. 



4-2.5 Summary 

The spectral, positional and temporal properties of the 
gamma-ray source 1FGL J1103.9— 5355 are inconsistent 
with originating from PSR J1103-5403. It seems most prob- 
able that 1FGL J1 103. 9-5355 is associa ted with PKS 1101- 
536 as described in iMahonv etafl I^OIOT ). Therefore we con- 
clude that the discovery of PSR J1103-5403 was entirely 
coincidental, despite the seemingly low probability. 



4.3 PSR J1604-44 

PSR J1604-44 has a spin period of 1.4 s, which is not typical 
of the pulsars detected by Fermi. Since its selection for this 
survey, 1FGL J1604. 7-4443 has been associat ed with the 
flat sp ectrum radio source PMN J1604-4441 (|Abdo et all 
boiObh . As discussed in Section [3TT1 the long integration time 
of our survey means that it is possible for us to detect one or 
more pulsars unrelated to the gamma-ray sources. Therefore 
we believe it is quite likely that this is a coincidental Galactic 
plane pulsar along the same line of sight as 1FGL J1604.7- 
4443 and PMN J1604-4441. Unfortunately the flux density 
of this pulsar is so low that the normal timing procedure is 
impractical with the Parkes radio telescope, and so we are 
unable to provide further details here. 



4.4 PSR J2241-5236 

4-4- 1 Radio timing 

The pulsar J2241-5236 has a spin period of 2.19 ms and a 
period derivative of 6.6 x 10 -21 ss _1 . This places the pul- 
sar firmly amongst the 'recycled pulsars', which have under- 
gone a period of accretion-induced spin-up from a compan- 
ion star. A full timing solution, and some standard derived 
parameters for J2241-5236 are given in Table [4] The en- 
ergy loss rate, E is 2.5 x 10 34 ergs -1 , well within the typical 
values for gamma-ray emitting pulsars. An oft-used gamma- 
ray detectability measure is E 1 ^ 2 jd 2 , which if we accept the 
DM-derived distance of 0.5 kpc, is 6 x 10 17 erg 1/2 pc~ 2 s~ 1/2 , 
typical of the known gamma-ray MSPs (|Abdo et alll2010ch . 

Figure [4] shows the polarisation profile of PSR J2241- 
5236, averaged over many hours of observation. The profile 
appears to be composed of several overlapping components. 
The pulse width at 50% of the peak is 68 /is and the pro- 
file shows some linear and circular polarisation, strongest 
towards the trailing edge of the pulse. 

The pulsar is in a circular orbit with an orbital period of 

3.5 hours. The mass function is 8.6 x 10 -7 Mq, which assum- 
ing a pulsar mass of 1.35 Mq, implies a minimum compan- 
ion mass of 0.012 Mq. Excluding planetary objects, this is 
the smallest mass function seen in pulsars outside of globular 
clusters. The standard explanation for the companion would 
be a white dwar f, from which mass has been stripped by 



be a white dwar t, from which m ass has been stripp ed by 
the pulsar wind dKluzniak et ail 1 19881 IPhinnev et~ailll988l ; 
Ivan den Heuvel fc van Paradiis Il98&t ). Similar to the PSR 
J0610-2100 i|Burgav et al.ll2006h . PSR J2241-5236 does not 
show any eclipse or DM variability throughout its orbit. It 
is likely therefore that the orbit is sufficiently inclined to the 
line of sight such that we do not 'see' the material removed 
from the companion. For this reason, we deem it unlikely 
that any Shapiro delay will be measured in this system. 

Given the high flux density and low DM of PSR J2241- 
5236 we have sought to determine the prospects for measur- 
ing proper motion and parallax effects, which would provide 
a model-independent distance to the pulsar. To do this we 
have simulated a 4 year timing campaign, observing every 
30 days, with 0.7^s RMS in the residuals. The real RMS 
in such a campaign will likely be dominated by unmodelled 
'timing noise', however this simulation does show the poten- 
tial for future timing of PSR J2241-5236. For our simulation 
we assumed the DM distance of 0.5 kpc and a velocity of 
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Table 4. Astrometric, rotational and orbital parameters from the 
timing of PSR J2241-5236. Values in parentheses are the nominal 
lcr uncertainties in the least significant digits quoted. 



PSR J2241-5236 



Right ascension, a 22:41:42.01850(5) 

Declination, <5 -52:36:36.2260(3) 

Dispersion measure, DM (cm~ 3 pc) 11.41085(3) 

Rotation measure, RM (rad m~ 2 ) 14(6) 

Pulse frequency, v (s" 1 ) 457.3101497559(3) 

First derivative of pulse frequency, v (s~ 2 ) — 1.388(10) x 10" 15 

Epoch of model (MJD) 55044.2 

Binary model BT^ 

Orbital period, P b (d) 0.1456722395(3) 

Epoch of periastron, T (MJD) 55044.1580909(6) 

Projected semi-major axis of orbit, x (lt-s) 0.02579537(14) 

Eccentricity <lx 10~ 5 

Validity of solution (MJD) 55044.5—55509.5 

RMS of residuals (/is) 0.7 

Clock correction procedure TT(TAI) 

Solar system ephemeris model DE405 

Coordinate time standard TCB 

Flux density at 1.4 GHz (S1400, mjy) 4.1(1) 

Width at 50% of pulse peak (Wso, /is) ... . 70 

Width at 10% of pulse peak (W w , /is) 130 

Distance* (kpc) 0.5 

Characteristic age (yr) 5 X 10 9 

Surface magnetic field strength, B (G) .... 1.2 X 10 s 

Rate of energy loss, E (ergs -1 ) 2.5 X 10 34 

Mass function (M Q ) 8.6844(4) x 10~ 7 

Minimum companion mass (Mq) 0.012 



assumes thelCordes fc Laziol J2002T) model of electron density. 
t |Blandford fc Teukolskvl dl976h . 

100 km s _1 . The simulated data were fit for frequency, fre- 
quency derivative, position, binary parameters, proper mo- 
tion and parallax. The error in the proper motion was less 
than 0.1 mas yr -1 , equivalent to a transverse velocity error 
of less than 1 km s~ . The parallax measurement gave an 
error of 0.4 mas, which equates to a distance error of 80 pc. 
Although real observations may well be adversely affected 
by additional sources of noise, this simple simulation indi- 
cates that PSR J2241-5236 has the potential to be studied 
in great detail. We intend to continue monitoring this pulsar 
for several years to measure the parameters such as paral- 
lax, as well as to investigate the pr ospects of using P SR 
J2241-5236 in a pulsar timing array (jriobbs et al.ll2010l ). 



4-4- 2 X-ray observations 

PSR J2241-5236 was covered by a programme of observing 
unassociated Fermi-LAT sources with the Chandra X-ray 
Observatory. On 2009 August 30, we observed this source 
for 20 ks with the ACIS-I instrument. The brightest X-ray 
source within the 1FGL error circle is at (J2000.0) R.A. = 
22:41:42.01, Decl. = -52:36:36.2 (90% radius = 0.3", sta- 
tistical only), less than 0.1 arcseconds from the position ob- 
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Figure 4. The 1369 MHz integrated pulse profile of PSR J2241- 
5236 clipped to a region round the pulse peak. Total intensity is 
shown as a thick line, with linear polarisation as a thin solid 
line and circular polarisation as a dotted line. Flux density is 
measured as a fraction of the pulse peak. The polarisation position 
angle variation over the pulse is shown in the upper panel. The 
polarisation position angles shown are not absolute. The inset 
figure shows the total intensity profile over the entire 360°. 

tained by radio timing. Therefore we are confident that PSR 
J2241-5236 is the source of the observed X-rays. 

We measured 86 counts (0.5 - 7 keV) arriving within 
5"of the pulsar position, whereas only 1.75 counts are ex- 
pected due to background. Therefore we model the spectrum 
using Cash statistics without subtracting the background. 
The X-ray spectrum was initially fit with an absorbed single 
black-body model with temperature, kT = 0.26 ± 0.04 keV 
(c-stat/d.o.f. = 174/254), but clearly left residual emission 
at <0.8 keV unmodelled. Adding a second black-body (kT = 
0.07 ± 0.04 keV) improved the fit (c-stat/d.o.f. = 167/252) 
and did not affect the temperature of the first black-body 
{kT = 0.31 ± 0.08 keV). The total 0.5 - 7 keV X-ray flux 
in the latter model was (3.1 jl5' 6 ) x 10 -14 erg cm -2 s _1 . 
The effect of Galact ic absorption (N n = 1.21 x 10 20 cm" 2 ; 
Kalberla et al.ll2005l ) in the fits was negligible, and the un- 
absorbed flux was 3.3 x 10~ 14 erg cm" 2 s" 1 . Assuming a 
distance of 0.5 kpc equates to a luminosity of ~1 x 10 30 erg 
s" 1 . This is less than 1% of the rotational energy loss rate, 
and so can be easily powered by the available spin down 
energy budget. 

4-4-3 Gamma-ray pulsations 

To search for gamma-ray pulsations, we have phase- 
binned the gamma-rays from 1FGL J2241. 9-5236 using the 
ephemeris obtained from the radio timing. To maximise the 
detection significance we searched over 20 values of radius 
and 20 values of low energy cut-off. With an extraction ra- 
dius of 1.2° and a low energy cut-off of E > 350 MeV we 
obtain a H-test parameter of 108, equivalent to a detection 




Figure 5. Phase aligned radio (solid line) and gamma-ray (his- 
togram with errors) integrated profiles of PSR J2241-5236. The 
gamma-rays from 1FGL J2241. 9-5236 were phase binned using 
the ephemeris obtained by radio timing. Two full periods are 
shown for clarity. 



significance above 8a (|de Jager fc B iisching 2010)). The best 
profile is shown in Figure[5] produced from the 1123 photons 
arriving between MJD 54684 and 55511 that matched our 
cuts. We are therefore confident that PSR J2241-5236 is the 
source of gamma-rays in 1FGL J2241. 9-5236. It is not clear 
if the gamma-ray profile has two or three components, how- 
ever we can measure the radio lag (i.e. the separation of the 
radio peak and the first gamma-ray peak) to be 0.15 ± 0.01 
in units of pulse phase. 



4-4-4 Gamma-ray spectral analysis 

We analysed the gamma-ray spectrum of 1FGL J2241.9- 
5236 using the same 21 months time interval from 2008 Au- 
gust 4 to 2010 May 4, as the 1FGL J1103. 9-5355 data set 
and we used the same analysis cuts and methods described 
in Section 14.2.31 We analysed a region of 10° radius around 
the radio position of the PSR J2241-5236, and modelled all 
of the 22 1FGL sources within 15° of the pulsar along with 
Galactic and isotropic diffuse emission as detailed in Sec- 
tion 14.2.31 There are no known pulsars in this region. Thus 
we did not use any cuts on the phase, and we modelled all 
sources by power law spectral shapes. The spectral param- 
eters for sources within the 10° of 1FGL J2241. 9-5236 were 
left free in the likelihood fit, and the spectral parameters for 
other sources were fixed to the values obtained from 1FGL 
catalogue. Due to the small number of photons per phase 
bin, we did not consider performing a phase-resolved spec- 
tral analysis. 

We modelled the spectrum of 1FGL J2241. 9-5236 with 
a power law with a simple exponential cutoff shape given by 



dN 
IE 



KE~ l exp 



\ ^cutoff / 



(2) 



where three parameters, the differential flux K, the pho- 
ton index V, and the cutoff energy -E cut off were allowed to 
be free in the fit. We found the best fit value of (1.62 ± 
O.llstat ± 0.03 sys ) for the photon index, a cutoff value of 




Energy (GeV) 

Figure 6. The gamma-ray spectrum of 1FGL J2241. 9-5236, the 
gamma-ray source associated with PSR J2241-5236. The butter- 
fly curve (solid lines) show the best-fit spectral shape of power law 
with exponential cutoff in the full energy range of 0.1 — 300 GeV. 
The flux points were obtained from independent fits in each en- 
ergy bin as explained in the text. 



sy S ) GeV, and a differential flux 
) xlO -12 MeV -1 cm-V 1 . The 



(4.0 ± 0.9 st at + 0.9sys - 0.6, 
of K=(7.7±0.7stat ±0.5 sys . 
integrated photon flux above 100 MeV was found as Fioo = 
(4.1 ± 0.4 stat ± O.lsys) xl0~ 8 cm~ 2 s~\ The errors quoted 
are the nominal ler statistical (first value) and systematic 
(second value) errors. The uncertainty in the LAT effective 
area is estimated to be 5% near 1 GeV, 10% below 0.1 GeV 
and 20% over 10 GeV. The resulting systematic errors on the 
three spectral parameters, propagated from the uncertain- 
ties on the LAT effective area, were calculated using a set of 
modified IRFs bracketing the nominal (Pass6 v3) one. The 
best fit spectrum in this full energy band is plotted as the 
solid curve in Figure 6. Only statistical errors were shown 
and taken into account in calculating the butterfly curve. 

The flux points in the figure were obtained from inde- 
pendent likelihood fits in each energy bin. For those fits, 
we assumed a power-law spectrum with index fixed at 2 for 
all the point sources and diffuse emission components. In 
this way, we obtained the flux values for all the sources in 
those energy bins, assuming that the widths of the energy 
bins were sufficiently small that all spectral shapes can be 
approximated with a power law in that range. The highest 
energy bin is determined by the highest energy photon from 
the source with a direction within the 68% containment ra- 
dius defined by the instrument point-spread function. For 
the energy bins where the source was not detected with TS 
> 9, we calculated an upper limit for the flux at 95% con- 
fidence level. In order to test the assumption of a cutoff 
shape in the spectrum of 1FGL J2241. 9-5236, we also fit- 
ted its spectrum with a power-law model. As expected for 
pulsar emission, this model was rejected by 6.8a using the 
likelihood ratio test (|Mattox et all) 19961 ) . 



5 CONCLUSION 

Observations of six gamma-ray pulsars discovered by the 
Fermi satellite were carried out using the Parkes radio tele- 
scope. We did not detect radio pulsations above a flux den- 
sity limit of ~ 0.04 mjy. Implications of these results are dis- 
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cussed in detail in lSaz Parkinson et al.l (|2010T ) and lRav et al.l 
j20ld ). 

In addition we have carried out deep radio observa- 
tions of eleven unassociated gamma-ray sources. Periodic- 
ity searches of these data have yielded the discovery of 
the millisecond pulsars J 1103-5403 and J2241-5236 and 
the long period pulsar J1604-44, and the detection of the 
MSP J1514-49. For the seven others, no radio pulsations 
were detected above a phase-averaged flux density limit of 
~ 0.04 mjy. This limit assumes a pulse duty cycle of 10%, 
no binary motion and that the pulsar is located at the centre 
of the beam. 

PSR J1604-44 is a weak, long period pulsar and we 
consider its association with 1FGL J1604. 7-4443 unlikely. 
Since its selection for observation, independent evidence has 
shown that 1FGL J1103. 9-5355 is as sociated with the fla t 
spectrum radio source PKS 1101-536 |Mahonv et al.ll201fj) . 
Furthermore, timing of the millisecond pulsar PSR J1103- 
5403 yields a position some 9' from the centroid of the 
gamma-ray source. Spectral and temporal analysis suggests 
that the gamma-ray source is indeed powered by an AGN. 
This demonstrates that with sufficiently deep observations 
it is possible to detect even rare objects by chance. 

PSR J2241-5236 has a spin period of 2.2 ms and orbits a 
very low mass companion with a 3.5 hour orbital period. The 
relatively high flux density and low DM of PSR J2241-5236 
makes it an excellent candidate for high precision timing ex- 
periments. When the photons from 1FGL J2241. 9-5236 are 
folded with the radio timing ephemeris, a complex gamma- 
ray profile is observed. The association of PSR J2241-5236 
with 1FGL J2241. 9-5236 is therefore clear. This is consistent 
with the gamma-ray spectral properties and the association 
with an X-ray source. 
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